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Abstract 

We study leptons in holographic composite Higgs models, namely in models possibly 
admitting a weakly coupled description in terms of five-dimensional (5D) theories. 
We introduce two scenarios leading to Majorana or Dirac neutrinos, based on the 
non-abelian discrete group S4 x Z3 which is responsible for nearly tri-bimaximal 
lepton mixing. The smallness of neutrino masses is naturally explained and nor- 
mal/inverted mass ordering can be accommodated. We analyze two specific 5D 
gauge-Higgs unification models in warped space as concrete examples of our frame- 
work. Both models pass the current bounds on Lepton Flavour Violation (LFV) 
processes. We pay special attention to the effect of so called boundary kinetic terms 
that are the dominant source of LFV. The model with Majorana neutrinos is com- 
patible with a Kaluza-Klein vector mass scale tukk ^ 3.5 TeV, which is roughly the 
lowest scale allowed by electroweak considerations. The model with Dirac neutri- 
nos, although not strongly constrained by LFV processes and data on lepton mixing, 
suffers from a too large deviation of the neutrino coupling to the Z boson from its 
Standard Model value, pushing ttikk ^ 10 TeV. 



1 Introduction 



The idea that the Standard Model (SM) Higgs might be a composite particle arising from a 
strongly coupled theory pQ has received considerable attention lately. One of the main reasons 
of this renewed interest comes from the observation that the composite Higgs paradigm is closely 
related to theories in extra dimensions [2] . This connection is particularly transparent in Randall- 
Sundrum (RS) models [3], thanks to the AdS/CFT duality |_[|. More precisely, certain theories in 
extra dimensions, including RS models, can be seen as a (relatively) weakly coupled description 
of a sub-set of 4D composite Higgs models. They consist of two sectors: an "elementary" 
sector, which includes the gauge and fermion fields of the SM, and a "composite" sector, which 
is strongly coupled and gives rise to the SM Higgs. The form of the couplings between these 
two sectors is not the most general one allowed by symmetry considerations only, but is more 
constrained. We denote in the following this more constrained class of models as Holographic 
Composite Higgs Models (HCHM). 

The flavour structure of HCHM has been studied in detail in the past mostly in the 5D 
context of RS models with fermion and gauge fields in the bulk and it has been shown to be 
particularly successful [5]. It automatically implements the idea of [6] to explain the hierarchy 
of the quark and charged lepton masses in terms of field localization in an extra dimension. 
Moreover, HCHM are equipped with a built-in GIM mechanism that goes under the name 
of RS-GIM [7] and automatically protects the SM fields from possibly large flavour violating 
interactions coming from the composite sector0 

Small neutrino masses and large lepton mixing are not easily accommodated in this set-up, 
because the large mixing potentially leads to excessive LFV. Neutrino oscillation experiments 
clearly show that the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix has a very 
peculiar structure well compatible with Tri-Bimaximal (TB) mixing [9]. There has been much 
progress in recent years in explaining TB lepton mixing and the absence of LFV interactions for 
charged leptons by means of discrete non-abelian symmetries. It is thus natural to apply such 
symmetries also in the context of HCHM in order to resolve the aforementioned problems. 

Aim of this paper is to introduce a class of HCHM where, thanks to a non-abelian discrete 
symmetry, lepton mixing is nearly TB, and at the same time bounds on LFV processes in 
the charged lepton sector are satisfied (see [10] for other proposals). The mass spectrum in 
the neutrino sector can be normally or inversely ordered. The pattern of flavour symmetry 
breaking is dictated by symmetry considerations only, without relying on extra assumptions 
[llj or specific mechanisms for the breaking of the flavour symmetry, such as the ones used 
in [T2J H3] (see also [H] ) in the case of A4 to reproduce TB mixing [TS] . We discuss the case 
of flavour symmetry breaking in the elementary and composite sectors to certain non-trivial 
subgroups of the original symmetry without advocating an explicit realization of the breaking^) 

1 Despite this protection mechanism, a CP violation problem is still present in the quark sector [3[8]. 
2 For 5D models, this is the flavour counterpart of breaking by boundary conditions, commonly used for gauge 
symmetry breaking. 
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In particular, no flavons or other specific sources of flavour breaking are present in our set-up. 
We consider in this paper the discrete group £4 x Z3. The group £4 has been shown [16] to be the 
minimal group giving rise to TB lepton mixing using symmetry principles only. The presence 
of an irreducible two-dimensional representation is another feature of S4. Such a representation 
allows to disentangle the symmetry properties of the third generation from the first two and is 
expected to be important when applying the flavour symmetry in the quark sector. 

We focus on two possible scenarios which only differ in the way SM neutrinos get a mass. In 
the first one, the SM neutrinos are Majorana fermions and the type I see-saw mechanism explains 
the smallness of their masses, with no need to introduce additional (intermediate) mass scales in 
the theory. In the second one, SM neutrinos are Dirac fermions and tiny Yukawa couplings are 
naturally explained by the ultra-composite nature of the right-handed (RH) neutrinos [17]. In 
both scenarios, the flavour symmetry is broken to Z2 x Z2 x Z3 in the elementary and to Zg^ in 
the composite sector. Note that the strength of this symmetry breaking is in general expected 
to be 0(1). In the composite sector for the charged leptons such a large breaking is actually 
favoured, because it allows to decrease the degree of compositeness of SM leptons, suppressing 
large deviations from the SM Ztt couplinglf] The breaking felt by neutrinos in the composite 
sector is instead required to be weak in the Majorana scenario, in order to not perturb too 
much TB lepton mixing. An alternative is to resort to an extra symmetry protecting neutrinos 
from being affected by the flavour symmetry breaking in the composite sector. On the contrary, 
flavour symmetry breaking in the composite sector can be large in Dirac models, provided that 
the tiny component of RH neutrinos in the elementary sector is flavour universal. 

After a general presentation of the basic 4D flavourful HCHM, we pass to construct two 
explicit realizations in terms of 5D warped models. For concreteness, we consider the HCHM 
where the Higgs is a pseudo-Goldstone boson, i.e. gauge-Higgs unification models [18]. The 
5D models are based on the minimal 5*0(5) x U(l)x gauge symmetry [19] , while the flavour 
symmetry group contains, in addition to the S4 x Z3 factor, model-dependent discrete abelian 
factors necessary to minimize the number of allowed (and often unwanted) terms. 

In the Majorana model, the leading source of flavour violation arises from so called fermion 
boundary kinetic terms (BKT), whose effect is analyzed in detail. The only sizable constraints 
come from lepton mixing, being LFV processes for charged leptons below the current bounds. 
We also argue that CP violating effects, such as the Electric Dipole Moments (EDM) for charged 
leptons, are negligibly small. Keeping the prediction of the solar mixing angle #12 within the 
experimentally allowed 3a range requires flavour symmetry breaking in the composite sector to 
be at most of 0(3% -j- 4%) for neutrinos, unless a Z2 exchange symmetry is present on the IR 
brane, in which case no constraint occurs. This Z2-invariant 5D model is surprisingly successful, 
simple and constrained, and essentially contains only one free real parameter and two Majorana 
phases! The model is compatible with the mass of the first Kaluza-Klein (KK) gauge resonances 

3 This is an important point, also because uncalculable contributions to LFV processes (and flavour preserving 
quantities as well) coming from higher dimensional operators are sub-leading with respect to the calculable ones 
only if the SM fields are sufficiently elementary. 
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being uikk ^ 3.5 TeV, which is roughly the lowest scale allowed by electroweak considerations 
(S parameter). The masses of all fermion KK resonances (charged and neutral) are always above 
the TeV scale. 

In the Dirac model the most significant constraint does not arise from LFV processes or 
lepton mixing, but from a too large deviation of the gauge coupling of neutrinos to the Z boson 
from its SM value, which is constrained by LEP I to be roughly at the per mille level. This 
bound is satisfied by taking ttlkk ^10 TeV, well above the LHC reach, with an 0(1%) tuning in 
the electroweak sector. The masses of charged fermion KK resonances are above the TeV scale, 
while in the neutral fermion sector potentially light (sub- TeV) states can appear. 

The structure of the paper is as follows. In section 2 we describe our set-up from a general 
effective 4D point of view both for the Majorana and Dirac models. In section 3 we briefly review 
the relevant operators entering in the LFV processes we focus on, radiative lepton decays l\ — > 
Z27, decays to three leptons l\ — > hhh an d /i — e conversion in nuclei. In section 4 we construct 
the 5D Majorana model, compute its mass spectrum in subsection 4.1, the deviations from 
gauge coupling universality in subsection 4.2, LFV processes and lepton mixing in subsection 
4.3 and estimate uncalculable effects in subsection 4.4. In section 5 a similar, but more concise, 
analysis is repeated for the Dirac model. We conclude in section 6. Three appendices are added. 
In appendix A basic definitions and properties of S4 are reviewed, in appendix B we report 
our conventions for the SO(5) generators and representations and in appendix C we write the 
detailed structure of the two form factors governing the charged lepton radiative decays. 

2 General Set-up 

We consider CHM with a non-abelian discrete flavour symmetry Gf = S4 x Z3. They consist of 
an "elementary" and a "composite" sector: 

£-tot — £-el ~\~ £-comp ~i~ £"mix • (2T) 

The symmetry Gf is broken in the elementary sector to Z2 x Z2 x Z3, where Z2 x Z2 C S4 
is generated by S and U, and in the composite sector to Zg , the diagonal subgroup of the 
external Z3 and Z3 C S4 generated by T (see appendix [A] for our notation and details on 
S4 group theory). We do not need to specify how the flavour symmetry breaking pattern is 
achieved. The term £ m ix governs the mixing between the two sectors and is assumed to be 
invariant under the whole flavour group Gf. This is our definition of HCHM in the following. 
We have two different classes of models, depending on whether neutrino masses are of Majorana 
or Dirac type. We will refer to the two cases as Majorana/Dirac models (or scenarios). 

2.1 Majorana Models 

The elementary sector is invariant under the SM gauge group and includes three generations 
of SM left-handed (LH) and RH leptons /£, and three RH neutrinos i/j^. Here and in the 
following Greek letters from the beginning of the alphabet denote generation indices; depending 
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on the context, a = e, /i, r or equivalently a = 1,2, 3. The LH leptons Z£ and the RH neutrinos 
u R transform as (3, 1) under S4 x Z3, while the RH leptons transform as (l,u; 2 ( a— where 
w = e 2ni/3 is the third root of unity. The elementary Lagrangian (up to dimension four terms) 
is taken to be 

c el = llibll + T%ibi a R + p%iBv% - \(y^M a ^ R + h.c.) , (2.2) 

where the superscript c denotes charge conjugation and M is the most general mass matrix 
invariant under Z2 x Z2 X Z3. In flavour space, it is of the form 

M = U TB M D U^ B , (2.3) 

with Utb the TB mixing matrix 



Utb 



1 *A 

3_ V 3 u 

l/l/l 
6 V 3 V 2 

,1 ,/I _. fl , 

\ V 6 V 3 V 2 / 



(2.4) 



and Mc a diagonal matrix. We use the notation A = "f^A^, for any vector A^. 

The composite sector is an unspecified strongly coupled theory, that gives rise, among other 
states, to a composite SM Higgs field. The latter may or may not be Goldstone fields coming from 
a spontaneously broken global symmetry. In absence of any interaction between the elementary 
and the composite sector, the SM fermions are massless. They gain masses, after ElectroWeak 
Symmetry Breaking (EWSB), by mixing with fermion operators \& belonging to the strongly 
coupled sector. The mixing Lagrangian £ m ix is 



where A is a high UV cut-off scale of the composite sector, ^f L , ^f R and ^?„ R are fermion 
operators of (quantum) dimensions 5/2 + 7^, 5/2 + 7^, 5/2 + 7^, transforming as (3,1), 
(ljW 2 ^ -1 )) and (3,1) under S4 x Z3, respectively. The mixing parameters A; L and X UR are 
flavour universal, while are flavour diagonal, but non-universal. For simplicity, we assume 
that all of them are real. Although strictly not necessary, we take jiljIiji ^ ^> so that these 
mixing couplings are irrelevant. To a good approximation, 11 and l R can be identified with the 
SM fields, with a small mixing with the strongly coupled sector. Integrating out the composite 
fermion operators and taking into account that C comp is invariant under Z 3 ^ only, gives the 
following charged lepton mass matrix (in left-right convention, ipiM^B.) 

- ^t^tKK) - Kvuk L Xt R 5 a ^Y R+llL , (2.6) 

where vh is the electroweak scale, fi is the O(TeV) scale at which the composite theory becomes 
strongly coupled and b a are 0(1) coefficients Q The hierarchy of the charged lepton masses 

4 The estimate (|2.6[) and the following are only valid for spontaneously broken CFT, which is the case mostly 
relevant for us. At a more qualitative level, however, our arguments apply to more generic holographic composite 
sectors, such as the ones in |20| . 
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naturally arises from the (///A) suppression factor in (|2.6p with a proper choice of anomalous 
dimensions 7^ and 7^. The coupling X UR is in general relevant and strongly mix with the 
composite sector. The latter gives the following contribution to the kinetic terms of v^. 

(2.7) 

with b a O(l) coefficients. When / j l/R < 0, the kinetic term in (|2.7p dominates over the 0(1) 
term (|2.2p present in the elementary Lagrangian, and it is more appropriate to say that z/j| are 
states in the composite sector with a small component in the elementary sector. When i/k are 
canonically normalized, the relevant coupling X UR in (|2.5p becomes effectively a constants The 
canonically normalized neutrino Dirac mass terms are of the form 



M D „ h L K R f^\~^R 1 f = p a . b Q v H Xi L fnyiL 

with b a 0(1) coefficients. Notice the crucial difference between the charged lepton (|2.6p and 
neutrino (j2.8|) masses. The former explicitly break the flavour symmetry, since I'jIr is not S4 x Z3 
invariant, while the latter do not, being I^vr an invariant. This implies that the coefficients b a 
vanish in the limit of exact S4 x Z3 symmetry, while b a = b, b a = b become flavour independent. 
Assuming a small breaking of the flavour symmetry in the neutrino sector, one can take b a ^ 6, 
b a « b and, independently of b m the Dirac neutrino mass terms ()2.8p become universal. We stress 
the importance of having a small breaking of the flavour symmetry in the neutrino sector but 
not necessarily for charged leptons, because the masses of the latter are already suppressed by 
their small degree of compositeness. Demanding a higher degree of compositeness, in particular 
for the r lepton, might result in too large deviations of its coupling to the Z from its SM value. 
Integrating out in this limit gives the following see-saw like neutrino mass matrix: 

M v>af) ~ w^^a) 2 ^'! 171 '^ 1 ^)^ ' (2 - 9) 

where we have again taken into account the scaling required to canonically normalize v^. Thanks 
to the latter, the factors b cancel from the final formula (|2.9p but, more importantly, we gain a 
crucial enhancement factor (/n/A) 27 "fl , without which the light neutrino masses would be far too 
small for Mb ~ O(Mpi), Mpi being the reduced Planck mass, considering that the Dirac mass 
terms are at most of 0{vjj)- No intermediate mass scale has then to be advocated for Mp. 

The mass matrix (|2.6p is diagonal in flavour space and no rotation of charged leptons is needed 
to go to the mass basis. On the other hand, the neutrino mass matrix (|2,9p is diagonalized by 
the matrix Utb (up to phases), which leads to the identification 

Upmns = Utb- (2-10) 

5 This can also be seen by solving an equation for the renormalization group flow of the couplings A |21l HQ] . 
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2.2 Dirac Models 



The elementary sector includes, like in the Majorana scenario, three generations of LH and RH 
leptons Z£, If^ in (3, 1) and (l,cj 2 ' a_1 * ) ) of S4 x Z3, respectively, and in addition we now have LH 
exotic neutrino singlets in (3, 1). The composite sector is assumed to contain two massless 
RH fermion bound states, singlets under Gsm, both in (3, 1) of S4 x Z3. One of them mixes 
with z>£ giving rise to vector-like massive neutrinos u a '. The remaining fermions, denoted by 
mix with some heavy vector-like states in the elementary sector. When the latter are integrated 
out, one is left with a tiny mixing mass term e between i/£ and z>£ in the elementary sector, 
which is flavour universalis The elementary Lagrangian (up to dimension four terms) is taken 
to be 

c el = itibit + i\ibi\ + i>tidi>t - e(^^ + h.c.) , (2.11) 

with M as in (|2.3p . The mixing Lagrangian C m i x is 



The operators , and are of dimensions 5/2 + 7/i, 5/2 + 7^, 5/2 + 7^^, transforming 
as (3, 1), (1, w 2 ( Q_1 )) and (3, 1) under S4 x Z3, respectively. The mixing parameters \ L and \p L 
are flavour universal, while are flavour diagonal, but non-universal. The charged lepton mass 
matrix is the same as (|2.6p . The operators ^f" R excite, among other states, the RH massless 
neutrino bound states that pair up with ££. The vector- like mass of v a depends on the nature 
of the coupling \„ L : 

m" ~ d a \y L ^(^-X 1 , for 7p L > 0, 



.A, 

m$ ~ d a Xo L n, for 7^ < 0, (2.13) 

where d a are 0(1) coefficients. When EWSB occurs, Yukawa couplings between \pp and 
induce mixing among v? and v a . When i> a are integrated out, one gets 

»2~%^yv", (2.i4) 



m?, VA 



V 



where d a are 0(1) coefficients. Plugging (|2.14p into the mass term in (|2.1ip gives the SM 
neutrino mass matrix 

M Vj(X p ~ e — M ap . (2.15) 

In the limit in which the masses and the mixing are universal, m" = m^, d a = d, d a = d, 
the mass matrix (|2.15p leads to TB mixing. The composite nature of the RH neutrino naturally 
explains the smallness of e and hence the actual SM neutrino masses |17j . 



6 We do not specify here how this mixing is achieved and how its flavour universality is guaranteed. We will 
see that the latter requirement can naturally be fulfilled in our 5D example, see subsection 5.1. 
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In both scenarios, the flavour symmetry Z3 , present in the composite sector, remains 
unbroken in the limit in which the neutrino mass term M in the elementary sector is proportional 
to the identity. Correspondingly, all tree-level flavour changing charged gauge boson interactions 
are vanishing in this limit. When M is not proportional to the identity, the latter are still 
negligible in the Majorana scenario, being suppressed by the masses of the heavy RH neutrinos, 
but can be sizable in the Dirac one, leading to processes such as fi — > ej. Tree-level flavour 
violating Higgs and neutral gauge boson interactions vanish in both scenarios. This summarizes 
the basic set-up of our Majorana and Dirac HCHM. There are of course several sub-leading 
effects that should consistently be analyzed. We have not performed such analysis, but have 
preferred to postpone their discussion to the explicit 5D models that will follow. We only 
comment here that a relevant source of flavour violation arises from the elementary sector, since 
the kinetic terms of the SM fermions (in the basis where £ m i x is Gf invariant) are constrained 
in general to be only Z2 x Z2 x Z3 invariant, rather than S4 x Z3 invariant: 



with Z\ = UtbZ^U^b, and Zf* a diagonal matrix. Similar considerations apply of course to 
and while the additional unbroken Z3 symmetry forbids flavour violating kinetic terms for 
Ift. As we will see, in 5D models the Z\ factors are mapped to BKT at the UV brane. 

3 Effective Field Theory for LFV Processes 

In this section we review, closely following [22] and their notation, the most relevant effective 
operators entering in LFV processes. The most experimentally constrained LFV observables 
are the radiative lepton decays l\ — > the decays to three leptons — > fohh an d the fi — e 
conversion in nuclei. Particularly relevant are the muon decays fi — > e"f and /i — > eee (/x — > 3e 
for short). These LFV processes are described by the following effective dimension 5 and 6 
operators 



Terms of the form {jji R eL){eLe R ) and (p-Le R )(e R eL), by a Fierz identity, are shown to contribute 
to $5 and g§, respectively. The first two terms contribute to fi — > e-y while all terms contribute 
to fj, — > 3e. One finds the following branching ratio for these processes: 



l L iDl L -> l L (l + Z{)iDl L 



(2.16) 



V2 



C e ff D m^ARfiRcr^eLF^ + m ti ,A L n L o llv exF^ + gi{^Re L ){e R e L ) + g2(^LeR){e L e R ) 
+ 5 , 3(AR7 /i eij)(eR7 At e jR ) + g<i(fiL^e L )(e L j^e L ) 

+ fi , 5(Afi7 /i ei?)(e L 7 At e L ) + 56(^L7 /i ei)(eR7 At e jR ) + h.c. . (3.1) 



4G F 



BR(fi -> &y) 



BR(/j, -> 3e) 




(3.2) 



8eRe[A R (2gl + g* 6 ) + A L (2g* 3 + g* 5 )] + 64e 2 ( 




y\A L \ 2 + \a r \ 2 ) 
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The fj,— e conversion in nuclei is more involved and described by an additional set of effective 
operators, that contain quark fields. The most relevant ones are the vector 4 fermion operators: 

~G^ CeffD 2 l(9LV(q)eLl^(iL + gRV(q)eRj' J 'fJ-R)q'y IJ ,q + h.c. . (3.3) 

q=u,d 

The branching ratio is given by (see [22] for more details) 

m 5 G 2 F 2 a 3 Z 4: 

B CO nv^N -> eN) ~ e// (l( 2Z + ^W^) + + 2JV)^ V(<0 | 2 

2 



+ | (2Z + N)g RV{u) + (Z + 2N)g RV(d) \ z j , (3.4) 

where Z and A" are the proton and neutron numbers of the nucleus, F p is the nuclear form 
factor, Z e ff is the effective atomic charge and T cap t is the total muon capture rate. A similar 
analysis applies to LFV processes involving the r lepton, see e.g. [23] for details. 

4 Explicit 5D Majorana Model 

It is useful to construct a specific 5D weakly coupled description of our Majorana scenario, in 
order to concretely address its phenomenological viability beyond possible estimates based on 
naive dimensional analysis only. We consider in the following a gauge-Higgs unification model 
in warped space \19\ 124] . As known, these models describe the sub-class of CHM where in the 
composite sector (a spontaneously broken CFT) a global symmetry G is spontaneously broken 
to a sub-group H, giving rise to a set of Goldstone fields including the SM Higgs field [25]. We 
consider the minimal symmetry breaking pattern SO (5) —> SO (4), leading only to the SM Higgs 
doublet. We use the conformally flat coordinates in which the 5D metric reads 

ds 2 = a 2 {z) (r] flu dx' i dx v - dz 2 ) = ( Vfll/ dx^dx u - dz 2 ) . (4.1) 

The UV and IR branes are located at z = R ~ 1/Mpi, where Mpi is the reduced Planck mass, 
and at z = R' ~ 1/TeV, respectively. The gauge symmetry in the bulk is 



and the flavour symmetry is 



Ggauge = 50(5) x U(l) x (4.2) 



Gflavour = 5 4 x Z 3 x Z' 3 X Z([ . (4.3) 



The gauge symmetry breaking is standard, with G gaugc broken at the UV and IR boundaries 
to Ggaug^uv = SU{2) L x U(l) Y and G gaugCi i R = 50(4) x U(l) x x P LR , where P LR is a LR 
Z2 symmetry, useful to suppress deviations of the couplings of fermions to the Z boson from 
their SM values [26]. The flavour symmetry is broken to Ofl avoU r,uv = Z2 x Z2 x Z3 x Z3 and 
avour, ir — Z^ ^ x Z3, respectively. In order to constrain the number of terms allowed at the 
UV and IR boundaries, two additional symmetries Z 3 and Z3 have been included. 
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Bulk 


UV 


IR 




Gflavour = £4 x Z3 X Z3 X Z3 


Gflavour,UV = Z 2 X Z 2 X Z3 X Z3 


Gfl a vour,IR = Z3 ' X Z' 3 




(3, l,w,w) 


(1,-1,1,0;) 
(-1,1,1,0;) 
(-1,-1,1,0;) 








(l.l,^- 1 )^) 






(3,l,w,l) 


(1,-1,1,1) 
(-1,1,1,1) 
(-1,-1,1,1) 





Table 1: Transformation properties of the 5D multiplets £j, £ e and £ u under Gfl avoU r and their 
decomposition properties under the subgroups Gfl avourj uv an d Gfl av our,iR m t ne Majorana model. 



The lepton particle content of the model consists of 5D bulk fermions only: one fundamental 
£i >a , one adjoint £ ejQ and one singlet representation £ v>a of 50(5), for each generation, all neutral 
under U(l)x (see [27] for a similar construction), 




/ x aL (+-) \ 
Vol(-\ — ) Z aL (-\ 

daL ( ) 

L>2,aL (H ) , L a L (H — 



(4.4) 



where the first and second entries in round brackets refer to the + (— ) Neumann (Dirichlet) 
boundary conditions (b.c.) at the UV and IR branes, respectively. We have written the 50(5) 
multiplets in (|4.4p in terms of their SU(2)l x SU(2)r decomposition, where [^1,^2] denotes 
the two components of the bi-doublet (2,2) with T^r = +1/2 (ipi) and T%r = —1/2 ("i/> 2 ). 
The SM LH lepton doublets arise from the zero modes of the 5D field L a i in the 5, the RH 
charged lepton singlets arise from the zero modes of e a R, T$r = — 1 component of the SU(2)r 
triplet in the 10, and the RH neutrinos arise from the zero modes of the singlet f^/jlll Notice 
that with the embedding (|4.4p . the LH SM charged leptons, originating from 5D fields with 
T3R = T%l = —1/2, are expected to have suppressed SM Z coupling deviations. In addition to 
the SM fields and their KK towers, the 5D fields (|4.4p give also rise to a set of exotic particles. 
In terms of SU(2)l x U(1)y, these are two doublets ii )Q .£ and L 2>Q £ with Y = 1/2, one doublet 
L a i with Y = — 1/2, two singlets u a L and D a i with Y = 0, one singlet x a L with Y = 1 and one 
triplet Z a i with Y = 0. The flavour properties of the fields (|4.4p are summarized in table [TJ 
Notice that the decomposition of the 3 of 54 into representations of the remnant group Z 2 x Z 2 
at the UV boundary implies a non-trivial basis transformation, see appendix 1X1 



7 The hypercharge Y and electric charge Q are given by Y — X + T;m and Q = T^l + Y. 
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The most general G gaug e ir x Gfl avoU r m invariant mass terms at the IR brane are 



'gauge 

1 



- Ar = (jj J ^2 ("4r, q ^Ll t aLL 2 ,aR + L aL L aR ) + mf R)Q i> a LV a R + h.c^j , (4.5) 

\ ' a=e,fj,,T 

all flavour diagonal. The only Gg aU g ej uv x G{[ aiVOUIj \jy invariant mass terms at the UV brane are 
Majorana mass terms for RH neutrinos: 

+ \ ™ UV , M + !fa + T^j (u eR + u^r + u tR ) 
+ ^ "iuv.r (v£ R - Kr) ~ v tr) + h.c. 

1 



2 VaR-Mvv,apv/m + h.c. (4.6) 



with 



Muv = U T Bmu V U^ B , (4.7) 

muv = diag (muv,e , m\JV,fi , m uv,r) and C/tb as in ()2.4p . Notice that the UV and IR localized 
mass terms are dimensionless. The phases of the IR mass terms mJ R a and mj^ a can be removed 
by properly re-defining the 5D SO (5) fields Q and ^ e>Q . We can also remove one of the three 
phases of the UV mass terms muv,a, so that, in total, the Majorana model contains just two 
phases. 



4.1 Mass Spectrum 

The mass spectrum of the theory (including all KK states) is efficiently computed using the 
so-called holographic approach [28] , which is also very useful to match the 5D theory to the 4D 
description given in section 2. As far as the lightest modes are concerned, however, simple and 
reliable formulas are more easily obtained using the more standard KK approach and the so 
called Zero Mode Approximation (ZMA), which we use in the following. The ZMA is defined 
as the approximation in which EWSB effects (i.e. Higgs insertions) are taken as perturbations 
and mixing with the KK states coming from Higgs insertions is neglected. The spectrum of 
the zero modes is then entirely fixed by the unperturbed zero mode wave functions and their 
overlap with the Higgs field. These unperturbed wave functions satisfy the new b.c. as given by 
the localized IR terms. As explained in [29], the localized UV Majorana mass terms, instead, 
must be considered as a perturbative mass insertion (like the Higgs) if one wants to recover a 
meaningful mass spectrum for the light SM neutrinos without taking into account mixing with 

8 Following a common use in the literature, we have omitted to write certain fermion terms, including terms in 
which the would-be fermion components with Dirichlet b.c. appear, because at the level of mass mixing the IR 
Lagrangian is only relevant in determining the modified b.c. (such as H4.8|l below) of the fields at the IR brane. 
One can detect the absence of such terms by noticing that the variation of the sum of the bulk and brane action 
at the IR brane does not vanish when the b.c. H4.8[) are imposed. 
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the KK states. Due to the wave function localization of zero and KK modes, as a rule of thumb, 
the lighter the zero mode masses are, the more accurate the ZMA is. 

Taking into account the localized IR mass terms (|4.5p . the IR b.c. for the non- vanishing 5D 
field components in ZMA are 



L a R = — TtiR, a L a R i L a i = m l m a L a L , z = R . 



m 



J 1 



R' 

z\ 2 / z 
R 



We get the following zero mode expansion: 

L aL (x,z) 
L aL (x,z) 
u aR (x,z) 
u aR (x,z) 
e a R(x,z) 



'Pa 



1R>a VR> Vi? 



i 



1 



-m 



IR,q" 



z \ 2 ( Z \ c l 



R> \RJ \R> 



1 



iR 7 \R 



R> 



f-c a e aR ( x ) > 



(4.8) 

(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 



where l^l(x), e^ R (x) and N^ R (x) are the canonically normalized LH lepton doublets, RH 
charged leptons and RH neutrino zero modes, respectively. We use the standard notation 



fc 



1 - 2c 



!-(#) 



l-2c 



1/2 



(4.14) 



where c = MR are the dimensionless bulk mass terms of the 5D fermions. We denote by q 
and c v the bulk mass terms of £; )Q and £ v ,a, constrained by the flavour symmetry to be flavour- 
independent. We denote by c a the remaining 3 bulk mass terms for £e,a* The parameters p a 
and <7„ are defined as 



Pa 



(t 1 



1 + \m 



IR,a\ 



f- 



f- 



(4.15) 



We take the unitary gauge for the 50(5) — > 50(4) symmetry breaking pattern in which the 
Higgs field wave function is (see appendix [B] for our 50(5) conventions) 



Al(x,z) = J^{ h -{x))=v H J^8^ =v H f H (z)6 & >\ 



with vh — 250 GeV. We find useful to introduce 



vh 



(4.16) 



(4.17) 
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where fn is the Higgs decay constant. It is defined as 

f B =*& = ? (418) 

g 5 R' gR'^\og{R'/R) 

In the second equality we have used the approximate tree-level matching between the 50(5) 5D 
coupling 55 and the SU (2) l 4D coupling g 



g 5 = g^/R\og(R'/R). (4.19) 

Flavour-independent bounds, essentially the S parameter in models with a custodial symmetry, 
constrain 1/R' > 1.5 TeV, corresponding to h < 1/3. 

By computing the wave-function overlap with the Higgs field, we get the following charged 
lepton mass matrix 

h 111} 

Ml,af3 = f Cl f-c a ^5 aj p . (4.20) 



'a 

As usual, the SM fermion masses are naturally obtained by taking c a < —1/2, in which case 
f- Ca are exponentially small and hierarchical. 
For the Dirac neutrino mass matrix we get 

The Majorana mass matrix in 4D is the one on the UV boundary. Taking into account the wave 
functions of the RH neutrinos, we get 

M^=(^r i it ^ (4,2, 

Integrating out the heavy Majorana fields N^(x), the factors a a cancel out and the actual form 
of the light neutrino mass matrix is, using (|4.7p . 



2R' Z 1 \ RJ Jp^ V m uv / a/3 Jpp 

For ~ TTiuv.a ~ C(l)) t ne size of the neutrino masses is mainly governed by the bulk mass 

term c„. The latter is essentially fixed to be 

c^w-0.36. (4.24) 

We have explicitly checked that the masses of the zero modes obtained in the ZMA (and treating 
the UV Majorana mass term as a perturbative mass insertion) are in excellent agreement with 
the exact tree-level spectrum. 

Let us consider the relation between the 5D model and the general 4D analysis performed in 
subsection [27TJ The strongly coupled sector is a CFT spontaneously broken at the scale fx ~ 1/R' 
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with a cut-off A ~ 1 /R. The anomalous dimensions appearing in (|2.5|) are uniquely fixed by the 
bulk masses of the 5D multiplets £ e and f„ [2T]PI 

m = |q + 1/2|-1, 7m = |Cc - 1/2| - 1, lvR = \c v - 1/2| - 1. (4.25) 

It is straightforward to show that, for q > 1/2 and c Q < —1/2, the (p/A) suppression factors 
appearing in (|2.6p arise from the factors f Cl and f- Ca defined in (|4.14p and that b a ~ ^ir q - 
With the value of cv taken as in (|4.24|) . the coupling X VR is relevant and /_ Cl/ ~ / Cl/ ~ 0(1)- 
The mass formula (|4.2ip is of the general form (|2.8|) . where b a ~ rrij^ a and b a ~ a a — 1. The 
latter factors, as expected, do not appear in the final mass formula (|2.9p . The p a are wave 
function normalization factors that take into account the contribution of the composite sector 
to the kinetic terms of the LH doublets, given by p a — 1. 

In the limit of an S4 invariant IR Lagrangian, mJ R a — > (so that p a — > 1) and a = mj^, 
the neutrino mass matrix (|4.23p leads to TB mixing. As we will see, bounds on gauge coupling 
deviations favour the region in parameter space where q is close to 1/2, in which case p a is equal 
to one (since f Ca > 1) to a reasonable approximation even for m[- Ra ~ 0(1). This accidental 
property allows us to also explore the region in parameter space where the flavour symmetry 
breaking in the charged lepton sector on the IR brane is large, while in the neutrino sector it 
remains small, namely m'( Ra = m^ R (l + Sm^^, with 6m'( Ra <C 1. 

Instead of assuming a small breaking in the neutrino sector and for the sake of reducing 
the number of parameters in the model, one might also advocate an accidental Z2 exchange 
symmetry present only in the IR localized Lagrangian, under which 

O a (x,R') ^ v a (x,R') . (4.26) 

If the symmetry ()4.26p is imposed, the IR mass parameters mj R a are constrained to be equal 
to ±1. Among the four inequivalent choices of ±1, we can take the universal choice m^ a = 1. 
Although not necessary, an analogous Z2 symmetry exchanging the two bi-doublets in the 5 
and the 10 of SO(5) (a single Z2 exchanging the bi-doublets and the singlets is also a viable 
possibility) might be advocated to also set wj Ra = 1. The resulting model can be seen as an 
ultra-minimal 5D model, with in total only 8 real parameters (5 bulk mass terms and 3 localized 
UV mass parameters) and two phases (contained in the UV mass parameters), 4 of which are 
essentially fixed by the SM charged leptons (c a , q), 1 by the overall neutrino mass scale (c„) 
and 2 by the neutrino mass square differences (two combinations of muv, Q ), leaving in this way 
just one free real parameter and two Majorana phases! We denote this constrained model by 
the "Z2-invariant" model. 

The mass spectrum of all the KK resonances is above the TeV scale. For instance, let us 
consider the Z2-invariant model and let us take h = 1/3, q = 0.52, c u = —0.365 as benchmark 
values. In this case, the lightest gauge boson KK resonances have ttikk ~ 3.5 TeV, the lightest 

9 The IR localized mass terms H4.5|l correspond to irrelevant deformations of the CFT and do not affect the 
anomalous dimensions computed in the limit of exact conformal symmetry. They however deform the mass 
spectrum of the CFT when finite cut-off effects are taken into account. 
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(negatively and positively) charged and neutral fermions have masses around 2 TeV, while the 
heavy Majorana neutrinos have masses around 10 13 GeV. 

4.2 Deviations from Gauge Coupling Universality 

In this subsection we compute the deviations from the SM values of the couplings of leptons to 
the Z and W bosons. In RS-like models such deviations can play an important role, since their 
expected order of magnitude for natural models with 1/R' > 1.5 TeV can be of the same order 
of magnitude or larger than the experimental bounds, which are at the per mille level. The size 
of the deviation is mainly fixed by the wave function profile in the fifth dimension of the 4D 
lepton. The more the field is UV peaked, the smaller the deviation is. On general grounds, one 
might expect sizable deviations for all the ZIlIl couplings and for the Ztrtr one. Deviations 
of the LH neutrino couplings Zu^ul should also be studied. The latter have indirectly been 
measured by LEP I and are constrained at the per mille level with an accuracy comparable to 
that for charged leptons, using the invisible decay width of the Z boson, under the assumption 
that this is entirely given by neutrinos [30j . 

An efficient way to compute these deviations, automatically summing over all the KK contri- 
butions, is provided by the holographic approach. In the latter, the effective 4D gauge couplings 
between fermions and gauge bosons are obtained by integrating over the internal dimension the 
5D gauge vertex, with the 5D fields replaced by bulk-to-boundary propagators and 4D fields. 
The main source of deviation arises from higher-order operators with Higgs insertions, which 
give a contribution of 0(h 2 ). Higher-order derivative operators are negligible, being suppressed 
by the fermion masses or Z boson mass and are OiMtR') 2 or 0{m z R') 2 , respectively. The 
momentum of all external fields (and hence of all bulk-to-boundary propagators) can be then 
reliably set to zero. In this limit, the computation greatly simplifies and compact analytic for- 
mulae can be derived. In the following we do not report all the details of our computation but 
only the final results. We define the 4D SM couplings gi^sM as 

9i,SM = Tl-Qsm 2 e w , (4.27) 

without additional factors of the coupling g or of the weak mixing angle 9\y- 

Let us first consider the LH charged leptons 11 . Given our embedding of /£ into 5D multiplets 
with T-jl = T^r, we simply have 

5gf L =gr L -g? L , S M = (4-28) 
and no deviations occur at all0 They occur for the RH charged leptons Z^. We get 



,(2 + 4q + (3 + 2c a )K RQ 



2-, 



2|m{ Ra | 2 (3 + 2c Q )(l + 2q) 



10 The coupling deviations above are denned in the field basis in which a completely localized UV fermion has 
SM gauge couplings, with no deviations. In this basis the fermion independent universal coupling deviation arising 
from gauge field mixing is encoded in the S parameter. 
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where M^ a are the charged lepton masses (|4.20p and it is understood that c a entering in (|4.29p 
are determined as a function of q, rw SRa and M^ a . Equation ()4,29p clearly shows that a small 
flavour symmetry breaking in the composite sector for charged leptons, i.e. <C 1, is 

disfavoured. Keeping q and M/ jQ , fixed, for small IR mass terms <5<7^ oc a \ 2 . It is 

intuitively clear that Sgf grow when the localized IR mass terms decrease, since one needs to 
delocalize more the RH leptons to get their correct masses, resulting in larger mixing with the 
KK spectrum and hence larger deviations. Alternatively, one has to decrease the value of q, 
increasing the degree of compositeness of the LH leptons. 

Let us now turn to the neutrino Z couplings. Since are embedded into SO (5) multiplets 
with T3L ^ non-trivial deviations are expected. In the limit R' 3> R and in the relevant 
range q > 0, c a < 0, c u < 0, we have 



h\l - 2c,) (4* - 2 + K Ri J 2 (2c " 3) + |mr R ,J 2 %^ 2 rj- 3) )^ 2c ' 



4(2q-3) (2c, 



1 + \ m ma\ 2 ( 2c i - ^))R'R 2ci + R{R'? Cl {l - 2c a ] 



(4.30) 



The couplings of the W boson to the LH doublets and their deviations from the SM values, 
denoted by 5g„ L i L i are computed in the same way. In the same limit as (|4.30|) . we find 



>l'L 



We demand that 



9l 



< 20/qo, 



Sg 



92 



< 40/00 



(4.31) 



(4.32) 



for LH and RH charged leptons and LH neutrinos. 

The LH neutrino deviations ()4.30p are mostly sensitive to q, requiring q > 0.49, with a mild 
dependence on the other parameters, while the RH charged lepton deviations are also very sen- 
sitive to the bi-doublet IR mass parameters, disfavouring small values of rn[ Ra . Independently 
of m| R a , we get an upper bound on q from the r lepton, q < 0.56. As in many warped models 
with bulk fermions, the region q ~ 1/2 is preferred by electroweak bounds. 



4.3 LFV Processes and BKT 

Due to our choice of discrete symmetries, no 5D operators that reduce to the operators appear- 
ing in (|3.ip and in (|3.3p are allowed in the bulk or on the IR brane. The flavour preserving 
dipole operators responsible for lepton EDM are also forbidden by gauge invariance. Operators 
associated with the couplings 54 and g§ in (|3.ip and gLV(q) m f|3.3j) are allowed on the UV brane, 
but their natural scale is 0{Mpf) and thus totally negligible. The operators in (|3.ip and in (|3.3p 
can only arise in the effective field theory below the KK scale, after the KK resonances have 
been integrated out. Their coefficients are then calculable. In absence of further corrections, 
tree-level flavour changing interactions among charged leptons mediated by neutral KK gauge 
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bosons and Higgs vanish, since all interactions and Yukawa couplings involving charged leptons 
are manifestly flavour diagonal (in contrast to what happens in generic RS models [31] ) : 

51-6 = 9LV(q) = 9RV(q) = ■ ( 4 - 33 ) 

Flavour violation occurs in the neutrino sector and hence radiative decays mediated by neutrinos 
and charged gauge bosons do not vanish, Al,Ar ^ 0. However, these are negligible, because 
effectively mediated only by heavy Majorana neutrinos. This is best seen by considering again 
the UV Majorana mass term as a mass insertion, but beyond ZMA, including all KK wave 
functions. The mass terms (|4.6|) can be written as follows: 

1 oo oo 

m=0 n=0 
(n) In) 

where f„aR are ^ ne KK wave functions for the fields and explicitly show that the Majorana 

mass matrix has rank 1 in the KK indices, for each flavour a. The heavy Majorana state N^ R 
is defined by the eigenvector in round brackets in ()4.34p . with non- vanishing Majorana mass. 
Suitable orthonormal combinations of define the "light" KK modes . The fields N^ R 
are not yet in their mass basis and (|4.34p is not diagonal in flavour space. However, since these 
fields are very heavy, we can integrate them out. In the limit of infinite mass, this implies 
setting N^ R = 0. Eventually, we see that the remaining terms in the Lagrangian involving the 
fields N^ft are flavour-diagonal with real coefficients. For finite mass, flavour and CP violating 
interactions are generated, but suppressed by the heavy Majorana mass and are completely 
negligible. 

It is important to study at this stage the impact of higher dimensional flavour violating 
operators in the model. These can only occur at the UV brane. The lowest dimensional operators 
of this form are fermion BKT. In principle all possible BKT allowed by the symmetries must 
be considered. In practice this is rather difficult to do, so we focus only on those BKT, whose 
presence with all others set to zero, causes flavour violation. From the table [H we see that Z3 
forbids the appearance of flavour violating BKT for £ e;a , while these are allowed for £j Q and 
£v a . There are in principle four possible flavour violating BKT at the UV brane, for Li jQ jj, 
VaRi v a R an d L a i. The KK expansion of fields with b.c. modified by both boundary mass and 
kinetic terms is quite involved. In order to simplify the analysis, we consider the BKT as a 
perturbation and treat them as insertions, like the Majorana mass terms. Namely, we take as 
b.c. for all fields the ones with vanishing BKT and then plug the resulting KK expansion into 
the BKT. This approximation is clearly valid for parametrically small BKT, but it is actually 
very good at the UV brane even for BKT of 0(1), as we will see (see [32] for an analysis of 
fermion BKT in warped models). Among the 4 BKT above, the UV BKT for L\ a R, v a R and 
u a fj are strongly suppressed (at least for the most relevant low KK modes), due to the form of 
the wave functions of these fields, and can be neglected. We are only left with 

Cbkt = L L (x, R)(RZi)iDL L (x, R) , (4.35) 
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where Z2 x Z2 constrains Zi to be of the form Z\ = JTrsdiag (z e i, Zni, z T i)U^ B . The coefficients 
z a i are dimensionless and their natural values are 0(1), although smaller values ~ l/(167r 2 ) can 
also be radiatively stable. If one assumes a small breaking of S4 — > Z2 x Z2 at the UV brane, the 
relative differences in the z a i can be taken parametrically smaller than ~ l/(16-7r 2 ). In presence 
of the flavour violating operators (|4.35|) . the couplings ()3.ip and (|3.3|) become non- vanishing. 

Let us first write down, in the mass basis, the relevant interaction terms of our 5D Lagrangian 
that give rise to the effective couplings present in (|3.1|) and (|3.3|) . We have 



Y, { C tJaLV- VlL + C? R l aR V-v iR + h.c.j + £ 9 g qV°q (4.36) 
i,a,V- q,V° 



+ Yl ( D t bI aLV°l bL + D a R % R V\ R ) + (YabUHlbR + Y b % R Hl bL ) 

a,b,V° a,b 

where a and b run over all charged leptons, q runs over the light SM quarks u and d, i runs 
over all the neutrinos, V~ and V run over all charged and neutral gauge fields, respectively. 
By "all" we here mean all species of particles, including their KK resonances. For simplicity of 
notation, we have omitted the implicit dependence of the couplings in (|4,36p on the gauge fields 
V~ and V°. The couplings in (|3.3|) depend on how the quark sector is realized in the theory. 
We assume here that up and down quarks are genuine 4D fields localized at the UV brane and 
singlets under the flavour symmetry. 

The coefficients Al and A R are radiatively generated and receive contributions from 3 dif- 
ferent classes of one- loop diagrams, 

Ar = 4T + 4 Z) + 4f } , Al = A L W) + A Z) + A™ , (4.37) 
where A^ L are the contributions due to the diagrams where a charged gauge boson and a 

(Z) 

neutrino are exchanged, A R L are the contributions due to the diagrams where a neutral gauge 

{ /■/ ) 

boson and a charged lepton are exchanged and A R j L are the contributions due to the diagrams 
where the Higgs and a charged lepton are exchanged. The explicit form of these coefficients is 
reported in appendix O 

The operators associated with the couplings g\-% are generated at tree-level by Higgs and 
neutral gauge boson exchange. By matching, we have 

222 2 

_ _ o m W y v * n — n m W V V _ _ \ ^ m W n ee _ \ ^ m W n ee 

9 1 ~ - 2 Z^T YeeY eii > 92 - , g 3 - ^ — —U R L) R , ff4-2^ -Z?T U L L > 

m H m H y m V° V m V° 

55 = £ ^DfD^ + ^Y ee Y e ^ , g 6 = £ ^D%Df + ^Y ee Y, e , (4.38) 

where mw and run are the masses of the SM W and Higgs bosons, respectively. The couplings 
in (|3.3|) are given by 

toVfc) = 4 £ ^ 9q D? , g RV{q) = 4 £ 4^ • (4.39) 

v o m V° yo m V° 



18 



Strictly speaking, the effective couplings appearing in (|3.1|) and (|3.3[) should be evaluated at the 
scale of the decaying charged lepton mass, while (|4.37p . ()4.38j) and (|4.39p give the couplings at the 
energy scale corresponding to the mass of the state that has been integrated out. Contrary to, 
say, non-leptonic quark decays, renormalization group effects in leptonic decays are sub-leading 
and can be neglected in first approximation. We can then directly identify the coefficients (|4,37D . 
(|4.38p and ()4.39p as the low-energy couplings relevant for the LFV processes. 

We have numerically computed the LFV processes by keeping, for each independent KK 
tower of states, the first heavy KK mode. For tree-level processes this approximation is quite 
accurate and should differ from the full result by 0(10%), as we have numerically checked by 
keeping more KK states. For radiative decays, the approximation is less accurate and might 
differ from the full result by 0(50%). This accuracy is enough for our purposes. If one demands 
a higher precision, a full 5D computation, as e.g. in [33], should be performed, although one 
should keep in mind that the limited range of validity of the effective field theory of 5D warped 
models puts a stringent bound on the accuracy one can in principle achieve. 

As we already said, all LFV processes are induced by the BKT (|4.35p . More precisely, LFV 
processes are induced by the relative differences in the z e i, Zui, z T i factors, since universal BKT 
simply amount to a trivial rescaling of the fields. Let us first give an estimate of the relative 
relevance of the couplings g\-g% and gL/RV{q)- They are induced by the tree- level exchange of 
Higgs and neutral gauge bosons, namely the SM boson Z and its first KK mode Z^ l \ the first 
KK mode of the photon the first KK mode of the neutral SO(5)/SO(4) fields A 3 ^ and 
A^ l \ the first KK mode of the 5D gauge field Z'W . The 5D fields Z, 7 and Z 1 are related as 
follows to the 50(5) x U(l) x fields W 3L , W 3R and X: 

95xW 3R + g 5 X , _ g^W 3R - g 5X X 

\Jai + gix ^al + glx 

cos OwW 3 l - sin 9 W B , 7 = cos 6 W B + sin 9wW 3 l , (4.40) 
with g§x the 5D coupling of the U(l)x field, determined in terms of 8\v- 

tan 2 9 W = / 5X 9 . (4.41) 
91 + 9tx 

Due to the IR-peaked profile of the KK wave functions, the leading effect of (|4.35p is to mix the 
LH zero mode fields among themselves. The main source of flavour violation clearly arises 
from LH fields. Since fermion Yukawa couplings are negligible, we have 

gi - 92 ^ 93 ^ 95 - gRV{q) - • (4.42) 

The LFV couplings in ()4.36p govern the size of the relevant effective couplings g±, g$ and 
9LV(q)- The dominant LFV effects arise from the rotation and rescaling of necessary to get 
canonically normalized kinetic terms. Before EWSB effects are considered, no flavour violation 
is expected from the SM Z boson by gauge invariance. The leading deviations arise from the 



B = 
Z = 
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gauge fields Z^ l \ and Z'^' . It is straightforward to derive a reasonable accurate estimate 
for the couplings D e ^: 

Df{Z'V) ~ -gCiZlU, (4.43) 

where <?z oc and gbulk are the BKT and bulk contributions to the gauge couplings, respectively. 
When EWSB effects are considered, LFV effects are transmitted to the SM Z boson as well. 
The resulting D e ^(Z) is suppressed by the mixing, but the latter is approximately compensated 
by the absence of the mass suppression factors appearing in the couplings gi (|4.38|) . Eventually, 
the SM Z boson contribution to LFV is of the same order of magnitude of that of the fields Z^ , 
7W and Z'^\ In (|4.43|) . Zi is the effective BKT felt by the zero mode, which is obtained by 
multiplying Z\ by the square of the zero-mode wave function (j4.9|) evaluated at the UV brane: 

( R ^ 1-2* n 1 /As 1 



(ZiU = - ^-=(^) Q/3 — = • (4.44) 
\R'J \JJ^ Jpp 



For c\ > 1/2, the factor entering in (|4.44|) becomes of 0(1), while it is exponentially small for 
ci < 1/2. For the relevant region where q ~ 1/2 and p a ~ 1, the effective BKT Zi is considerably 
smaller than Z\. For q = 1/2 + 5, at linear order in 6, we have 

Zi * ( log" 1 E+S^Z-Q^ + S^Zi. (4.45) 

The effect of the BKT on the LFV is naturally suppressed. This is the main reason why most of 
the parameter space of our model successfully passes the bounds imposed by LFV processes. The 
suppression factor (|4.45|) also explains why the approximation of treating the BKT as insertions 
is valid even for 0(1) BKT at the UV brane. 

Let us now consider the couplings Al and ^4^. It is immediately clear from the more com- 
posite nature of the muon with respect to the electron that Al <C Ah, so in first approximation 
Al can be neglected. Higgs and neutral gauge boson mediated contributions A R and A R are 
also negligible, and the dominant contribution A^ R ^ arises from the charged gauge bosons with 
Neumann b.c. at the UV brane, namely the SM W boson and its first KK excitation wj^ . It 
turns out to be rather difficult to derive an accurate analytic formula for A^^ since neutrino, 
charged lepton and gauge boson Yukawa couplings significantly contribute to the branching ra- 
tio. An order of magnitude estimate can be obtained by focusing on a definite contribution 
that is always one of the dominant ones, although not the only one. It arises from the Yukawa 
couplings between the SM neutrinos an d the RH singlet fields N l aR , the combination of 
N^ R and N^ R orthonormal to the heavy Majorana fermions N^ R . It is relevant because these 
Yukawas are sizable and N l a is typically the lightest fermion resonance in the model. We get 
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Figure 1 : Branching ratio of /i — > e~f and /i— e conversion in Ti as a function of the UV BKT Sz = z^i — z e i . 
The continuous (red) and dashed (black) lines in the left panel are the actual [33] and the expected future 
bound [35] given by the MEG experiment. The (red) line in the right panel is the current bound as given 
by the SINDRUM II experiment [55]. The plots refer to the Majorana model with TOj Rq = 1, cj = 0.52, 
c v = —0.365, h = 1/3 and normal neutrino mass hierarchy. The IR masses m\ R a arc random numbers 
chosen between 0.05 and 1.5 for m l IR e „ and 0.5 and 1.5 for m\ R T (blue points) or all set to one (green 
diamonds). The masses totjv.q are chosen such that the lightest neutrino mass is mo = 0.01 eV and the 
best fit values [37] of the solar and atmospheric mass square differences Am s 2 | = 7.59 x 10 -5 cV 2 and 
Am 2 tm = 2.40 x 10~ 3 eV 2 are reproduced using (|4.23p . corrected for the effect of the BKT. 

where c is an order 1 coefficient and Y is the approximate flavour universal value of the Yuakwa 
coupling in the original basis of fields, before the redefinitions needed to get canonically normal- 
ized kinetic terms. We plot in figure Q] the bounds arising from fj, — )■ ej and /i — e conversion 
in Ti (the most constraining case) as a function of Sz = z^i — z e \ = 3(Zi) efl . Both processes 
depend quadratically on 5z, as expected from (|4.43|) and (|4.46p . As can be seen from figure 
[U the IR masses mj R do not play an important role, provided that m,j R is large enough, as 
required by 5gJ . Thanks again to the suppression factor appearing in (|4.45p . the branching 
ratio is almost always below the current limit of 2.4 x 1CP 12 for \Sz\ < 1. Using however the 
future bound expected from the MEG experiment of 10~ 13 , we find that \5z\ is constrained to 
be less than 0.25. The decay to three leptons fj, — > 3e and radiative r decays are always well 
below the experimental bounds and are not reported0 We also performed an analysis for larger 
mo and for both, normal and inverted, neutrino mass orderings, with results identical to those 
shown in figure 1. 

Let us finally consider the bounds arising from lepton mixing, assuming vanishing phases. 
As we have already mentioned, in order to avoid too large deviations from TB lepton mixing, 
the IR localized neutrino mass terms iTi^ a should be taken close to universal. Parametrizing 
m ma i n ^ ne following way: m,j Re = mj R , rn^ = m^il + Sm^), mJ Rr = mj^f^l we can 

1 Notice that, due to the smallness of the couplings of the SM leptons to the neutral KK gauge bosons for 
c; > 1/2, the contribution of A), to (|3.2p is comparable to that given by the couplings 174 and g%. 

12 We have chosen this particular parametrization of m^ R a , since in this way all mixing angles are subject to a 
deviation linear in Sm\ n and neither accidental cancellation nor accidental enhancement of the coefficient of the 
linear perturbation is encountered. 



21 



analyze neutrino masses and mixing arising from the light neutrino mass matrix in (|4.23p in 
an expansion in 5m^. We neglect the effects of BKT in the following and take q = 0.52 and 
h = 1/3 so that the parameters p a are universal to a good approximation. For normally ordered 
light neutrinos with itlq = 0.01 eV and mass square differences given by the experimental best 
fit values [37], the mixing angles turn out to be 

sin #i3 0.05 l^mj^l , 

sin 2 #23 ~ 2 + °- 82 6m m > 

sin 2 0i2 « - - 1.58 <5mf R , (4.47) 
3 

showing that the requirement of having sin 2 #12 in the experimentally allowed 3cr range [37j . 
0.27 < sin 2 612 < 0.38, leads to the constraint 

- 0.03 <Sm^< 0.04. (4.48) 

At the same time sin 2 #23 remains within its la range, 0.475 < sin 2 #23 ^ 0.533. The reactor 
mixing angle sin 2 #13 takes as maximal value 4 x 10 -6 , well below the current and prospective 
future bounds. These statements are in agreement with our numerical results. The validity 
of the expansion in Sm^ strongly depends on tjiq. For instance, by taking Sm^ = 0.1, the 
perturbative expansion in Sm^ breaks down for uiq > 0.03 eV. 

We also performed a study for inverted mass hierarchy. In this case the above perturba- 
tive expansion is not valid for any value of ttiq. From the numerical results we see that large 
corrections to the solar mixing angle always arise, whereas the atmospheric mixing angle gets 
small corrections and the reactor mixing angle remains always very small. The large devia- 
tions of #12 can be easily understood by noticing that for inverted neutrino mass ordering the 
relative splitting between the two heavier light neutrinos is in general small compared to the 
scale y/ m\ + Am 2 tm > 0.049 eV. Thus, the angle #12 associated with the mixing in this almost 
degenerate sub-sector is subject to large deviations from its initial TB value even for very small 
deviations tfrnj^ from universality. As a consequence, the latter have to be as small as possible in 
the case of inverted neutrino mass ordering, which is most naturally achieved in the Z2-invariant 
model. 

In summary, in the case of normally ordered light neutrinos and a rather small mass scale mo, 
deviations from universality of m^ Ra are admissible up to the level l^m^pj < 0.04. Generically, 
the solar mixing angle, which is the most precisely measured one up to date in neutrino oscillation 
experiments, turns out to be the most sensitive one to corrections. For a neutrino mass spectrum 
with inverted hierarchy, the most natural situation is the one in which an additional accidental 
Z2 exchange symmetry on the IR brane renders the mass terms mj R a universal. The deviations 
of 612 and #23 are well under control in the Z2-invariant model for all values of mo and both types 
of neutrino mass hierarchy (with sin 2 #23 in the experimentally allowed la range and sin 2 612 in 
the 2a range). The angle #13 is in this case always constrained to be very small, sin 2 #13 < 10 -6 , 
and cannot be detected. 
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4.4 Uncalculable Corrections and r Decays 

Contrary to the operators in (|3.ip . where only two flavours appear, LFV operators involving 
simultaneously three different flavours are not constrained effectively by our choice of discrete 
symmetries. Dimension 8, 4 fermion 54 x Z3 (Zg^) invariant bulk (IR localized) operators 
reducing to flavour violating dimension 6 LL, RR and LR/RL operators can be constructed. 
Among these, the ones of the form (fTfi)(eTiJ,), (er//)(eTr), with T = 7 M ,7 5 , and their hermitian 
conjugates, can directly mediate the r decays r — > e2[i and r — > fj2e. The branching ratio for 
these decays is of order 10 -8 [38]. The size of the couplings of these operators, uncalculable 
within the 5D theory, can be estimated by using naive dimensional analysis. For IR brane 
operators (bulk operators give roughly the same result) we get 

Kuv9lR 3 (§) 4 m(z = X) c K Uv g\R'f log f[ f Cn , (4.49) 

n=l 

where c n = —c a for RH leptons, c n = c\ for LH leptons, Kfjy is an 0(1) dimensionless coupling 
and in the second equality we have plugged in the zero mode wave function of the SM leptons 
;(°) for the 5D fermion fields ^. The most stringent bounds arise from the LL operators, since 
/_ Ce ,/_ Cfl < f Cl . By matching with ([3J]), we get 

/R'\ 4 

g e Jf ~ 2k uv {m w R'f log (-) J] f Cn . (4.50) 

n=l 

We demand that g$} < 2 x 10" 4 @ For k uv ~ O(l), 1/R' > 1.5 TeV and 1/R ~ M PU we get 
the following lower bound for q: 

Q>0.5. (4.51) 

Notice that flavour preserving dimension 8 operators are also potentially dangerous, contributing, 
e.g., to the deviation from the SM values of the couplings of leptons to the vector bosons. From 
a quick estimate, we find that the bound (|4.5ip is more constraining. Summarizing, demanding 
that the uncalculable contributions coming from higher dimensional operators are sufficiently 
suppressed results in a bound on the degree of compositeness of the SM leptons. 

5 Explicit 5D Dirac Model 

In this section we provide an explicit 5D gauge-Higgs unification warped model realizing the 
Dirac scenario outlined in subsection 12.21 The model is very closely related to the Majorana 
model of section so we focus on the key differences between the two. The gauge symmetry 
and its breaking pattern is the same as before, while the flavour symmetry is slightly different: 

^flavour = Si X Z 3 X Z 5 X Z3 , (5.1) 

13 Notice that in r decays, the formula (|3.2p gets a suppression factor ~ 0.18 that takes into account of the 
hadronic contribution to the total decay width. This explains the factor 2 in 2 x 1CP 4 . 



23 





Bulk 


UV 


IR 




Gflavour = 54 X Z3 X Z5 X Zg 


Gflavour.UV = Z 2 X Z 2 X Z3 X Z3 


GflavourJR = Z^ X Z5 




(3, 1,0/5,0;) 


(1,-1,1,0;) 
(-1,1,1,0/) 
(-1,-1,1,0/) 










(u/ 2 (^),o/ 5 ) 




(3,l,o/f,o/) 


(1,-1,1,0/) 

(-1,1,1,0/) 
(-1,-1,1,0/) 


(o/ 2 ^- 1 )^ 2 ) 



Table 2: Transformation properties of the 5D multiplets £ e and under Gfl av0 ur and their 
decomposition properties under the subgroups Gfl avour] uv an d Ofl avour] iR in the Dirac model, uj^ 
is the fifth root of unity 0/5 = e 2m / 5 . 



broken to Gfl avourj uv = Z 2 x Z 2 x Z3 x Z' 3 and Gfl avouri iR, = Z 3 x Z5 at the UV and IR branes, 
respectively. Like in the Majorana model, Z5 and Z 3 are included to constrain the number of 
terms allowed at the UV and IR boundaries. 

The particle content and b.c. for the fields are identical to those in the Majorana model, 
with the only exception of a crucial flip in the b.c. for the singlet neutrino v in the 5: 

v aL (-+) -> VaL (+-) • (5.2) 

The flavour properties of the fields are summarized in table [2j Notice that the discrete symme- 
tries forbid the appearance of any bulk or boundary Majorana mass term. 
The invariant mass terms at the IR and UV branes are 

/ R \ 4 3 

- £lR = ( "^7 J ^2 m IR,a \Ll,aLL2,aR + ^qL^qR^) + h.C. 
\ ' a=e,fi,r 

-£UV = V aL M\JV,a/3V/3R + h.C. (5.3) 

with A^tjv as m (|4.7p . The phases of the IR masses m l IRa can still be absorbed by re-defining 
the 5D 50(5) fields £ ejQ and one of the three phases contained in A4uv through re-phasing the 
fields £v t a- Again, we are left with two non-trivial phases coming from the UV mass terms. 



5.1 Mass Spectrum 

The KK expansion in the ZMA of the doublets L a L and e a R is identical to (|4.9p . (|4.10|) 
and (|4.13|) and gives rise to the same charged lepton mass matrix (|4.2U|) . 

The KK expansion of neutrinos is of course different. The IR b.c. are not affected by mass 
terms, while the UV b.c. read 

V a R = Mw,al3^l3R, V aL = -Mjjy^UpL (5-4) 
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O^'TWyii) \«) KS) f-^sU^N^r.), (5.5) 



and lead to the following canonically normalized zero mode expansion 
u aR (x, z) = -= (-|) (^j) " f- Cv {UTB) a p-=Nfl{x) , 



K/3 

where 

By computing the wave function overlap with the Higgs field, we get the neutrino mass matrix: 
= ^fcf-c [ W ) -^{UtbU^-— , (5.7) 



where p a is defined as in (|4.15p . Thanks to the factor (R/R') in (|5,7p . the correct order of 
magnitude for neutrino masses is naturally obtained by choosing 

Cl/ -cj«0.8. (5.8) 

The only source of deviation from TB mixing in (|5,7p is given by the factor p a , which should be 
contrasted with the situation in the Majorana model, where the deviations are given by p a and 
the neutrino mass terms rn^ Ra . 

Let us consider the relation between the 5D model and the general 4D analysis performed 
in subsection 12.21 The anomalous dimensions of tyf R , *$>f R L are the same as in the Majorana 
model. Since ul arid ul belong to the same 5D bulk multiplet, we have "Jcl — Iil- The states 
denoted by v a in (|2.13p are the lightest KK vector-like states of the tower of modes coming from 
v a and u a in the 5D model. Their masses are determined as the zeros of a certain combination 
of Bessel functions and are approximately flavour independent. For q > 0.44, we have 

«R') 2 - , 1 " 4C ' | |2n - . r . (5.9) 

(n i l\ (R'\ 2c l- l l-t i |muV, Q! | 2 (l-2Ci) ^ V ' 

V C « 2> \R) \J~ T{cu-l/2)\l-2c v ) ) 

The coefficients d a appearing in ([2.13P are correspondingly flavour independent in first approx- 
imation. Along the lines of [21], the parameter e defined in subsection 12.21 can be seen to arise 
from the mixing between two heavy elementary fermions and *&r of opposite chiralities with 
a RH massless bound state vr of the CFT, all in 3 of S4. Omitting flavour indices, the relevant 
Lagrangian is 

C = $>{id - A)tf + -£—(* L i/ a + vr^l) + (hM* R + h.c.) (5.10) 

with ^y u R = \c u + l/2\ — 1 the anomalous dimension of vr and c a universal 0(1) coefficient. If 
we assume that A is flavour independent, all the non-trivial flavour dependence is in the mass 
term M, which is of the form (|2.3p . Integrating out the fields gives at leading order 

~ ; ° ,, vr. (5.11) 
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Rescaling vr — > pP vR+ vr to effectively get canonical kinetic terms for a free fermion field, gives 



The form of the coefficients d a introduced in (|2,14|) will be determined in subsection 15.21 We 
anticipate here that they are flavour independent, coming from £4 invariant bulk interactions. 
It is important to notice that for 7^ = 7^ = q — 1/2 (taking q + 1/2 > 0), (|2.14p shows a 
non-decoupling effect. For 7^ > the explicit (p/A) suppression factor in (|2.14p cancels the 
one coming from m^. For 7^ < 0, the mass of u is unsuppressed, but the LH leptons are mostly 
composite and one has to perform a field rescaling to get the canonically normalized LH field i>l, 
as in (|2.7p . Its effect is again to compensate for the explicit factor (/i/A) in (|2.14p . For any yp^, 
then, we get vl ~ 0(h)vL, a result that leads to unsuppressed deviations of neutrino couplings 
to the W and Z from their SM values, as shown below. The non-trivial factor K a in (|5.7p comes 
in the 4D picture from corrections to the kinetic term of vr we have neglected, appearing when 
the heavy fields \I/ are integrated out. They are completely negligible, given the suppression 
factor appearing in (|5.6p . The factors p a , as in the Majorana model, encode corrections to the 
kinetic term of coming from the composite sector. Summarizing, the mass formula (|5.7p is a 
particular realization of the more general expression (|2.15|) where all deviations from TB mixing 
are naturally suppressed. 

In the limit of an S4 invariant IR Lagrangian, m[ R — > 0, the neutrino mass matrix (j5.7[) 
leads to exact TB mixing. However, the factors p a disfavour composite LH leptons, because the 
more these states are composite, the smaller mj Ra should be to keep p a ~ 1. Bounds on gauge 
coupling deviations favour the region in parameter space where cj < 1/2. Given that TB lepton 
mixing and (|4.51|) favour q > 1/2, the region q ~ 1/2 is again the one of interest. 

The mass spectrum of the neutral KK resonances in the Dirac model differs from that in 
the Majorana model mostly for the presence of the light states v. For the benchmark values 
h = 1/3, Cj = 0.52, (|5.9p gives m- > 200 GeV. The masses of the next-to-lightest charged and 
neutral KK fermion resonances (taking m{ RQ between 1/2 and 3/2) are slightly below 2 TeV, 
so approximately comparable to the spectrum found in the Majorana model. The KK gauge 
boson masses are obviously identical in the two cases. 

5.2 Deviations from Gauge Coupling Universality 

The realization of the SM charged leptons in the 5D Dirac and Majorana models is identical, 
so (|4.28p and (|4.29p continue to apply. In the Dirac model, u a L are still embedded into 50(5) 
multiplets with T^l 7^ T$r, so non-trivial deviations are expected. The holographic analysis 
is complicated by the presence of the 4D singlet fields i) a L(x,z = R) that should be kept and 
eventually integrated out0 Omitting intermediate steps, one simply gets 



Recall that in the holographic approach the 5D fields evaluated at the UV brane can directly be identified 
with the fields in the elementary sector defined in section 2. 




(5.12) 




(5.13) 
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independently of any parameter. As anticipated below (|5.12p . a non-decoupling occurs in 
At leading order in h, (|5,13p leads to the following universal deviation 





h 2 




1 




9S L 


2 ' 


9 »lIl 


~~ 2 





(5.14) 



By demanding (|4.32p . we get the universal bound 

h<— (5.15) 
~ 10 

independently of q. In light of the bound (|5.15p . the Dirac model appears to be fine-tuned at 
0(1%) level, unless one advocates exotic hidden physics that is responsible for a fraction of the 
invisible partial width of the Z boson. 



5.3 LFV Processes and BKT 

Several considerations made in the Majorana model continue to apply in the Dirac case. The 
form of the interaction Lagrangian is the same as in (|4.36p . and the matching given by (|4.38p 
and $OgD still holds. The analysis in ()4T39]) - (f4T4"5]) is valid also here0 The bound on q (|4"3Tp 
coming from UV uncalculable corrections also applies here. 

In contrast to the Majorana model, radiative decays mediated by neutrinos and charged 
gauge bosons are no longer negligible, even in the absence of BKTj^J Interestingly enough, in 
this case we have been able to find a reasonable analytic formula for see (|5.18p . working 

in the flavour basis where Yukawa couplings are treated as perturbative mass insertions. Given 
the difficulty of finding such formulae, we report in the following some details on how (|5.18|) 
has been derived. We still adopt a KK approach and keep, for each 5D fermion field, only the 
first KK resonance. Even in this approximation, an analytic computation is complicated by the 
large number of fields that are present. The most important point to note is that the flavour 
violation comes from the singlet fields and arising from the expansion of the 5D fields 
(|5.5p . The zero modes due to their ultra- localization towards the IR brane, are effectively 

decoupled and can be neglected. Among the massive KK gauge bosons, the leading contribution 
comes from the charged gauge field in the SO(5)/ SO(4) coset, since it directly couples 
Na^ to the SM leptons. We can then safely neglect the SU(2)l x SU(2)r massive gauge fields 
WjP and Wft\ We neglect the tiny neutrino Yukawa couplings, their only effect being the 
rotation of the SM neutrinos with Utb- Charged lepton and gauge Yukawa mixing are in first 
approximation negligible. The relevant Lagrangian terms are the following: 

C D N^Y Nv U TB lT ] + 0) WU TB lf 0) + NP&AWlfV + h.c. (5.16) 

v2 



A similar non-decoupling effect has recently been noted in 5D models in flat space, see (3.15) of |20] , 
16 Notice that in the Dirac model, in principle, we might have LFV BKT on the IR brane. They arise from the 
RH KK neutrinos that, in analogy to the zero mode fields (|5.5[) . contain Utb in their expansion. However, this 
effect is indirect, and driven by the UV mass terms muv,a- Their impact on the model is sub- leading. We have 
numerically checked it in the BKT insertion approximation. 

17 The same is valid for CP violating effects, that we have not studied in the Dirac model. 
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(a) 




(b) 




Figure 2: Leading one- loop graphs contributing to fj, — > e*y in the flavour basis for (a) W and 
(b) AW exchange. The crosses represent Yukawa coupling insertions. 



where the flavour index has been omitted. The Yukawa couplings Y^ u are flavour non-diagonal 
with roughly the following structure: Y^ u ~ U^ b Yq^ u + 5Y^ U , where Yq : n u is a number and 
8Yn v a matrix in flavour space with |<5Y/v;,| <C |Yo,a^^| • The gauge couplings gj^ are also flavour 
violating and have the approximate form g^ ~ ^tb9q + &9 L ■> where g$ is a number and 8g L a 
matrix in flavour space with \ Sg L \ <C \g$ |. It turns out that the leading contribution to Ar comes 
from the first term in square brackets in A^\ see (jC.lj) . Indeed, the potential enhancement 
of the second term coming from the muon mass in the denominator is compensated by the 
smallness of the Yukawa coupling responsible for a non- vanishing RH coupling Cf R . The leading 
contributions coming from the W and A^ exchange are depicted in figure [2j Notice that no 
Yukawa insertion in the loop is needed in the diagram (b), because the relevant gauge interactions 
are already flavour violating. The computation of the two diagrams gives: 

A w ~ ie Y o,Nu Sm NM a aW -5ieg^5g L m 2 w 

R nc_2 _„3 ' R or. /o 2 n 2 — 2 ' V • / 



96vr 2 m 3 M (1) ' H ' 32^/3- 



taken Sg L all equal in flavour space. Terms proportional to m N (i) 5Yn u are sub-leading and have 
been neglected in A^ . In A^ W we keep the leading terms in the expansion m N (i) <C rn A (i 



where 6m N (i) = m (i) — m (i) is the mass splitting before EWSB and for simplicity we have 

(!)■ 

Indeed, in most of the parameter space, due to the chosen b.c. for the singlet in the 5, the 
neutrinos Na (which should be identified with the fields v a defined in subsection 2.2) are 
sensibly lighter than the 50(5)/50(4) gauge field A^\ In particular, for q > 1/2, m N (i) become 
very light. Roughly speaking, it turns out that \5g L \ < 1 5m N (i) /m N (i) | and |Yo,iVi/| m w, so 
that Aft^ /A^ 7 ~ ( m NW / m Aw) 2 ^ 1 an d the dominant contribution comes from the exchange 
of the SM W boson. Expanding (|5.9[) up to 0(|muv,a| 2 ) and using the ZMA formula (|5.7p . we 
get the following estimate for the branching ratio: 

Am s 2 ol log" 2 (?) (Y QyNu \i(R'\ic v -i 



BR{Li -> erf) ~ ^ 



6(1 + 2 Cl )T(c v + 3/2)/| [\ R> J \RJ 



(5.18) 



The Yukawa coupling Yq^ u depends of course on the input parameters as well, but there seems 
to be no simple expression for it. We have checked, by comparison with the full numerical 
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Figure 3: Branching ratio of fi — > ey and fi—e conversion in Ti as a function of the UV BKT Sz = z^ — Zei 
in the Dirac model. The continuous (red) and dashed (black) lines in the left panel represent the current 
[33] and the expected future bound [33] given by the MEG experiment. The (red) line in the right panel is 
the experimental bound as given by SINDRUM II [36] . The plots refer to the Dirac model with q = 0.52, 
c v = 1.33, h = 1/3 and normal neutrino mass hierarchy. The IR masses m[ R a are random numbers 



chosen between 0.05 and 1.5 for m m e and 0.5 and 1.5 for mj R T (blue points) or all set to one (green 
diamonds). The masses totjv.q are chosen such that the lightest neutrino mass is mo = 0.01 eV and the 
best fit values [37] of the solar and atmospheric mass square differences Am s 2 | = 7.59 x 10 -5 cV 2 and 



2.40 x 10 3 eV 2 are reproduced using (|5.7[) . corrected for the effect of the BKT. 



computation, that (|5.18p is accurate at the O(10%) level. The branching ratio crucially depends 
on the values of c v and q . 

When the UV BKT in (|4.35p are considered, the branching ratio of f/, — > ey receives extra 
contributions of the form (|4.46p that for \5z\ > 0.02 dominate. Unfortunately, it is not simple to 
derive a reasonably accurate analytic expression for BR(fi — > ey) in this case. We plot in figure [3] 
the branching ratio of /U — > ey and fi — e conversion in Ti for the Dirac model for c\ = 0.52 and 
c v = 1.33, setting all phases to zero. As shown in subsection I5.2[ the deviation of g" L from its 
SM value puts a strong bound on h, h < 1/10, see (|5.15p . We take here a value of h = 1/3 in 
order to compare the results for the LFV processes in the Majorana model with those in the 
Dirac model^£| The bounds are mainly governed by the UV BKT, with a very mild dependence 
on m[ R a . As can be seen from figure El both processes depend quadratically on 5z, as expected 
from (|4.43p and (|4.46|) . and the relative difference \5z\ of the UV localized BKT is constrained 
to be smaller than 0.15 in order to pass the actual MEG bound of 2.4 x 10" 12 . This becomes 
smaller than < 0.05 to pass the expected future MEG bound BR(/j, — > ej) < 10~ 13 . For such 
small values of 5z, cancellations between the contribution (|5.18p and the one associated with 
the UV localized BKT can occur and further suppress BR(fi — > e^/), see figure [3] The results for 
BR(fi -> 3e) and B com; (mTi — y eTi) are automatically below the current experimental bounds, 
as soon as BR(fi — > &y) is below the new limit set by the MEG Collaboration. Radiative r 

18 As we mentioned, one might take the optimistic point of view that the bound (|5.15[) might be a signal of new 
exotic hidden physics rather than a true bound for the model. Anyhow, for h < 1/10, no bound would arise from 
BR(fi —t ey) or _B con „(/iTi —t eTi), since both scale as h 4 , and the points in figure [3] should be scaled down by 
two orders of magnitude. 
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decays, r — > ^7 and r — > e~f, have branching ratios < 10~ 9 . The branching ratio of /i — > ey in 
the Majorana model is about two orders of magnitude smaller than the one in the Dirac model, 
for equal values of c\. In contrast, B conv {jj!Y\ — > eTi) is similar, being governed in both cases by 
the same tree-level FCNC. 

Concerning the lepton mixing angles, we find sin 2 #23 well within the experimentally allowed 
la range and sin 2 #12 still within the 2a range |37| . The value of sin 2 #13 is smaller than 10~ 8 . As 
already discussed in detail in the case of the Majorana model with non-universal masses m'( Ra , 
an inverted neutrino mass hierarchy is disfavoured, because the solar mixing angle receives 
in general too large corrections, while the atmospheric mixing angle still remains within the 
experimentally allowed la range and sin 2 #13 < 10~ 8 . A noteworthy effect in the Dirac scenario 
is the rather large deviation of the lepton mixing matrix Upmns from unitarity. For example, 
in the case of the set of parameters used to generate the plots shown in figure El we checked 
that the diagonal elements of U pmns Upmns an d UpmnsU pmns can deviate up to 0.05 from 
one. Their off-diagonal elements are in general much smaller. The non- unitarity of Upmns 
is associated with the non-decoupling of the light states Na\ This is in sharp contrast with 
the results found in the Majorana model in which the deviation from unitarity of Upmns 1S m 
general less than 10 -3 for the diagonal elements of U pmns Upmns and UpmnsU pmns - 



6 Conclusions 

We have introduced a class of 4D HCHM based on the non-abelian flavour group S4. x Z3, 
where lepton masses can be naturally reproduced and nearly TB lepton mixing is predicted. 
Both Majorana and Dirac neutrinos can be accommodated. A small breaking of the flavour 
symmetry for charged leptons is disfavoured in the composite sector, typically leading to a too 
large deviation of the coupling of the r to the Z from its SM value. The latter observation 
is linked to the choice of representations of the discrete flavour group used for the LH and 
RH charged leptons and needed to forbid large flavour violating effects. It applies to more 
general constructions based on different flavour groups. The breaking of the flavour symmetry 
for neutrinos in the composite sector can be large in the Dirac model, whereas it must be small 
in the Majorana model to suppress too large deviations from TB mixing. 

We have also constructed two explicit realizations of our framework in terms of 5D gauge- 
Higgs unification theories. We have computed in detail the relevant bounds coming from LFV 
processes in the charged lepton sector and shown that no significant constraints arise in both 
models. In the Majorana model, all the spectrum of fermion resonances is above the TeV scaler^ 9 ! 
while in the Dirac case light (sub-TeV) neutral fermions appear and are responsible for a too 
large deviation of the coupling of neutrinos to the Z from its SM value. A particularly economic 
and successful Majorana model can be constructed by postulating a Z2 exchange symmetry on 
the IR boundary which protects neutrinos from being affected by the flavour symmetry breaking. 
In both models, Majorana and Dirac, two CP phases are present. We have not studied in detail 

19 This should be contrasted with [13], where light fermion resonances appear. 
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their effects on the lepton EDM, but have argued that in the Majorana model these are expected 
to be negligibly small. 

Overall, the 5D Dirac model performs worse than the 5D Majorana model but, of course, 
this does not necessarily imply that more natural HCHM based on the Dirac scenario cannot 
be constructed, rather we might have missed to find a better representative. 

Note Added 

During the final stages of the preparation of this paper new experimental results have been 
released by the T2K [39| and MINOS [40] Collaborations indicating that #13 = is disfavoured 
at the level of 2.5a and 89% confidence level, respectively. Subsequently, three groups [4TJI42JI13] 
have performed a global fit of the available neutrino data finding at different levels of significance 
#13 ^ 0. The strongest indication of #13 7^ is found by [41] at a level of (more than) 3a, while 
the analysis in [43] shows that the mixing angle #13 is still compatible with zero at the latter 
level. The best fit value of sin 2 #13 is 0.01 -^0.02 in all three analyses |41 [ I42 |, . If such sizable 
value of sin 2 #13 will be confirmed in the future, our models (with Dirac and Majorana neutrinos, 
respectively) become disfavoured, because they generically foresee small values of sin 2 6*13 below 
10 -4 without additional (new) sources giving rise to #13 7^ 0. However, this does not rule out 
HCHM with flavour symmetries (broken in the manner as proposed by us) in general, since 
other mixing patterns, see e.g. [44] . apart from TB mixing, can be implemented as well |45| . 
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A Group Theory of 64 

In this appendix we briefly recapitulate the group structure of S4 and discuss the decomposition 
of the S4 x Z3 representations under the subgroups Z 3 ^ (preserved in the composite sector /on the 
IR brane) and Z2 x Z2 x Z3 (preserved in the elementary sector/on the UV brane), respectively. 
S4 is the permutation group of four distinct objects and is isomorphic to the symmetry group O 
of a regular octahedron. It has 24 distinct elements and five real irreducible representations: 1, 
1', 2, 3 and 3', out of which only the two triplets are faithful. We define S4 with the help of three 
generators S T and lE which are of the following form for the five different representations: 

20 S4 can also be defined in terms of only two generators. 
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1:5 = 
1' : 5 = 

2 : S = 
3: 5 = 

3' : 5 = 

and fulfill the relations 
S 2 




T 
T 




U 
U 

U 

u 
u 



1 . 



\0 co 




1 

\oio t 



1 , T d 



1, u 2 



1 



(ST) 3 = 1 , {SUf = 1 , {TUf = 1 , {STU) A = 1 . 



(A.l) 



Note that S and T alone generate the group A4, and similarly, that the two generators T and U 
alone generate the group S3. The character table of S4 can be found in, e.g., [36]. The Kronecker 
products are of the form 

1 x n = /j, V At , l'xl' = l, l'x2 = 2, 

l' x 3 = 3' , 1' x 3' = 3 , 

2x2 = 1 + 1' + 2, 2x3 = 2x3' = 3 + 3', 

3x3 = 3' x 3' = 1 + 2 + 3 + 3', 3x3' = l' + 2 + 3 + 3'. (A.2) 

The Clebsch Gordan coefficients can be found in, e.g., [U]. For (01,02,03)*, (01,02,03)* ~ 3, the 
invariant under S4 is of the form C1C1+C2C3+C3C2. Note that the choice of T being complex for the 
real representations 2, 3 and 3' leads, for (</>i,</>2)* ~ 2, (tp\, ip^, ipa)* ~ 3 and (V4; ^' 2 , V^)* ~ 3', 
to conjugate fields transforming as (</>2,<fii) t ~ 2, (tpl , ~ 3 and (i/i'*,^,^)* ~ 3 ' 
The decomposition of the S4 x Z3 representations under zi is given by 



(1, 


uP) - 


■+ w J 




(!' 


w*) - 


■» u; j 




(2, 


<J) - 


■» ^ 




(3, 


iJ) _ 


■» W j 


+ ^+2 + £ ^+l 


(3' 


w*) - 


■» U; j 


+ ^'+ 2 +^' +1 



(A.3) 



for j = 0, 1, 2. Since the generator T is diagonal in the group basis chosen by us, we can easily 
see that, e.g. for ~ (3,1), tpi transforms as 1 under Zg D \ ip2 as uj 2 and ^3 as w . The 
decomposition of the S4 representations under the subgroup Z2 x Z2 generated by S and U, 
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with SU = US, isO 

1 -> (1,1) 
1' -> (1,-1) 

2 -> (1,1) + (1,-1) 

3 -> (1,-1) + (-!,!) + (-!,-!) 

3' -> (1,1) + (-!,!) + (-!,-!) (A.4) 

where (±1, ±1) indicate the transformation properties under the two Z2 factors of Z2 x Z2. 
Since 5 and U are not diagonal, the decomposition of the S4 representations under Z2 x Z2 is 
non-trivial. For fa ~ 2 we get 

-^(&±fo)~(l,±l), (A.5) 

for a triplet ^ ~ 3 

-^(^l+^2 + ^ 3 )~(l,-l), ~= (V>2 - $}) ~ (-1, 1) , -^(2^i -^2-^3) ~ (-1,-1) 

(A.6) 

and similarly for ^ ~ 3' under S*4 

-L(^ + ^ 2 + ^)~(l,l), * ( 2 ^-^ 2 -^)~ (-1,1), (V 2 - ^3) ~(-l,-l)- 

(A.7) 

B 50(5) Generators and Representations 

We list here the explicit choice of 50(5) generators and SU(2)l x SU(2)r embedding used in 
the paper. Denoting by 

% = = (B.l) 

the 10 anti-symmetric generators of 50(5), where a, b = 1, ... ,5 label the generators and i,j 
their matrix components, we take 

*i = -^ 23 +* 14 ), tl = - l -{t^+f% ti = - l -{^+t'% 

t a = --^=t a5 , a = 1,2,3,4. (B.2) 

In this basis, t]] 2 ' 3 generate SU{2) L , t^ 2 ' 3 generate SU(2) R and t 1 ^ 1 G 50(5)/50(4). A 
fermion multiplet ^5 in the 5 of 50(5) decomposes as 5 = (2,2)0(1, 1) under SU(2) L x SU(2) R 



1 The external Z3 factor remains unbroken in this case and is omitted in the following. 
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and can be written as follows: 



1 

V2 



where 



q± 



f i{u + — d-) ^ 
-{u + + 
i(u- + d+) 
u_ — d + 
\ \/2n J 

u±\ 



(B.3) 



d±. 



(B.4) 



are the two doublets with T3R = ±1/2, respectively, forming the bi-doublet, and n is the singlet. 
A fermion multiplet *i in the 10 of 50(5) decomposes as 10 = (2, 2) © (1, 3) © (3, 1) under 
SU(2)l x SU(2)r and can be written as follows: 



V2 Y 
t 1 + it 2 



V2 



f 1 - it 2 



V2 



V2 



-u + + 



t 1 - it 2 



V2 



d- + 



t 3 + it 4 
V2 



-u- + 



t 3 - it 4 



V2 



-d+ , 



(B.5) 



where =t , (fP form the SU (2) ^ triplet , x^iX the SU(2)r triplet, and u±, <i± are the components 
of the bi-doublet, as defined in ()B.4|) . 



C Expressions for Al and A 



R 



We report here the explicit expressions for A^rJl an d -^#/£> i n terms of the couplings 

defined in (|4.36p in the main text: 



A 



(Z) 
R 



64ir 2 * — ' m;,. 
i,v- V 

-ie \ - m? 



64vr 2 



.4 



-ze \ - mr 



647r 2 * — ' m'zr 

a ti 



M, 



Til tj V ''i//. ' 



with = MHm\_, zf = M^/m^o, z% = M 2 /m 2 H , m w and m# the SM and Higgs 
masses, and 



1 



10 - 43z + 78z 2 - 49z d + 18z d log z + Az~ 



6(z- l) 4 



(C.2) 
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fz(z) = g^^u ( 8 - 38z + 39z2 ~ 18z2 lQ g z ~ Uz3 + 5 ^ 4 ) » 



6(z 

(T-l)8 



= 7 l 1 „ (4-3z + 6zlogz-z 3 ) , (C.3) 



/ff(*) = _\ )4 (2 + 3z + 6z log z - 6z 2 + z 3 
gff (*) = (z "^3 (3 + 2 log z - \z + z 2 ) . (C.4) 
All the expressions above are in agreement with the results of [48] . 
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